Reduced and projected two-particle entanglement at finite temperatures 
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We present a theory for two-particle entanglement production and detection in mesoscopic con- 
ductors at finite temperature. The entanglement of the density matrix projected out of the emitted 
many-body state differs from the entanglement of the reduced density matrix, detectable by current 
correlation measurements. Under general conditions reduced entanglement constitutes a witness for 
projected entanglement. Applied to the recent experiment [Neder et al, Nature 448 333 (2007)] on 
a fermionic Hanbury Brown Twiss two-particle interferometer we find that despite an appreciable 
entanglement production in the experiment, the detectable entanglement is close to zero. 
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The last decade has witnessed an increasing interest 
in generation and detection of entanglement in meso- 
scopic conductors 0, Q . Entanglement is an ubiquitous 
quantum effect, it describes correlations between parti- 
cles that can not be accounted for classically. A better 
understanding of entanglement of elementary charge car- 
riers, or quasiparticles, is therefore of fundamental inter- 
est. Due to controllable system properties and coherent 
transport conditions, mesoscopic conductors constitute 
ideal systems for the investigation of quasiparticle en- 
tanglement. In a longer time perspective, the prospect 
of quantum information processing using spin or orbital 
quantum states of individual quasiparticles provides ad- 
ditional motivation for such an investigation. 

To date quasiparticle entanglement has remained ex- 
perimentally elusive. However, recently an important 
step was taken towards a demonstration of entanglement 
in mesoscopic conductors. Based on the theoretical pro- 
posal [3| for a fermionic two-particle interferometer (2PI), 
see Fig. [TJ Neder et al 0| were able to demonstrate inter- 
ference between two electrons emitted from independent 
sources. In perfect agreement with theory, the interfer- 
ence pattern was visible in the current correlations but 
not in the average current. Under conditions of zero de- 
phasing and temperature, the part of the emitted state 
with one electron in each detection region A,B would be 
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where 1,2 denote the sources. The wave function \^> s ) is 
maximally entangled, it is a singlet in orbital, or pseudo 
spin, space. However, in the experiment [4| a reduced 
amplitude (~ 25%) of the current correlation oscillations 
was observed, suggesting an important effect of both de- 
phasing and finite temperature. This raises two interest- 
ing and interrelated questions: are the electrons reaching 
the detectors at A and B entangled and if so, can this 
two-particle entanglement be unambiguously detected by 
measurements of currents and current correlators? 

In this work we provide an affirmative answer to both 
these questions. We present a general theory for two- 



particle entanglement generation and detection in meso- 
scopic conductors at finite temperatures and apply it to 
the 2PI. Under very general conditions a nonzero entan- 
glement of the reduced orbital density matrix, accessible 
by current and current correlation measurements [H, is 
shown to be a signature of finite entanglement of the den- 
sity matrix projected out from the emitted many-body 
state. In other words, finite reduced entanglement con- 
stitutes a witness for finite projected entanglement. For 
the 2PI-experiment, while the projected density matrix 
is found to be clearly entangled, the reduced, observable 
density matrix is only marginally entangled. 
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FIG. 1: (color online) a) Schematic of the two-particle in- 
terferometer (2PI) [j| with beam splitters C,D and biased 
(grounded) contacts 1,2 (3,4). Detector regions A and B 
(yellow shaded) contain beam splitters and grounded con- 
tacts ±. b) Amplitude-temperature plot for the 2PI with 
semi-transparent splitters C,D. Above/to the left of the plot- 
ted lines the entanglement of the projected (reduced) density 
matrix is finite, C p > (C r > 0), and a Bell inequality is 
violated, Wgf > 2 (STg? x > 2). The red ring denotes the pa- 
rameters of the experiment 0|, showing that while the emitted 
state is clearly entangled, C v > 0, it is barely detectable by 
current and current correlation measurements, C r sa 0. 



We are interested in finite temperature orbital [6| en- 
tanglement in a general mesoscopic system, shown in 
Fig. [3] A theory for entanglement production in non- 
interacting pj mesoscopic conductors at finite tempera- 
ture was presented by Beenakker [l|. At a given energy, 
only the component of the emitted many-body state with 
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one particle in detector region A and one in B has nonzero 
entanglement. Formally the entanglement of the emitted 
state, here called projected entanglement, is quantified in 
terms of the two-particle density matrix projected out 
from the many-body state. In an experiment, while this 
two-particle density matrix in principle can be projected 
out by local operations/measurements and classical com- 
munication between A and B, it can not be directly ac- 
cessed by standard measurements of currents and current 
correlators The reason for this is twofold: 

First, at nonzero temperatures it is not only the bi- 
ased source reservoirs which emit particles but also the 
grounded source reservoirs and the detectors do. As a 
consequence, there is a finite amplitude for emitted states 
with two-particles at A and/or at B. These unentangled 
states contribute to currents and current correlators, re- 
sulting in a detectable, effective state with suppressed 
entanglement. Second, the current and current correla- 
tors provide information on the energy integrated prop- 
erties of the many-body state but not on the emitted 
state at each energy. This lack of energy-resolved infor- 
mation leads to a further suppression of the detectable 
entanglement. Clearly, these consequences of the ther- 
mally excited Fermi sea constitute generic problems in 
mesoscopic conductors. 

As a remedy for these finite temperature read-out 
problems it was suggested to work with detectors at very 
low temperatures [J. Recently, Hannes and Titov Q in- 
vestigated entanglement detection at finite temperatures 
via a Bell Inequality violation. In order to overcome the 
problem with detectors emitting particles they proposed 
to introduce energy filters, such as quantum dots at the 
drains. However, these schemes [3, Q would lead to ad- 
ditional experimental complications in systems which al- 
ready are experimentally very challenging. In this work 
we take a different route and investigate what informa- 
tion about the projected entanglement can actually be 
deduced from current and current correlation measure- 
ments. It is known Q that such measurements allow for 
a complete, tomographic reconstruction of the reduced 
orbital two-particle density matrix. At the focus of our 
investigation will thus be the relation between the pro- 
jected entanglement and the entanglement of the reduced 
density matrix, called the reduced entanglement. 

To provide a physically compelling picture, we first in- 
vestigate entanglement generation and detection in the 
2PI. Thereafter a formal derivation for a general meso- 
scopic system is presented. We consider the 2PI shown in 
Fig. [I] with source reservoirs 1 and 2 biased at eV while 3 
and 4 as well as the detector reservoirs are grounded 
All reservoirs are kept at the same temperature T. The 
reflcctionless source beam-splitters C and D have trans- 
parencies Tc — 1 — Rc and To = 1 — Rd respectively. 
Note that all electrons impinging on the detectors are 
emitted by the source reservoirs; scattering between the 
detectors is prohibited in the 2PI-geometry. 



In the 2PI-experiment, working with semi-transparent 
splitters Tc = Tjy = 1/2, a two-particle Aharonov-Bohm 
(AB) effect Q was observed in the current cross corre- 
lations Sa+b+- For finite temperature and dephasing, 
theory |lOj predicts 

Sa+b+ = -e 3 V/(4h)H [l- 7 sin & ot ] (2) 

where H — cothx — X/x with x = eV/2kT, 7 is a phe- 
nomenological decoherence parameter ranging from 1 for 
a fully coherent system to for an incoherent one and 
4>tot, up to a constant, is the AB- phase. The applied bias 
7.8/j.V and the estimated temperature lOmK in the ex- 
periment yield H = 0.78. A direct comparison to Eq. 
© then gives the oscillation amplitude H'y = 0.25, i.e. 
7 = 0.32, a substantial dephasing. To determine the two- 
particle entanglement of the emitted state we first calcu- 
late the projected (unnormalized) density matrix p p (E) 
at energy E. Using the formal similarity of the 2PI and 
the reflection- less, non spin- mixing splitter studied in [l[, 
Eqs. (B9) - (B13), we get 

p p (E) = (1 - ffft [ XP f aa + (1 - x)¥'1 (3) 

where x = e ~ 2x an d / = 1/(1 + e E / kT ) and /y = 
1/(1 + X eE ^ kT ) the Fermi distribution functions of the 
grounded and biased source reservoirs respectively. We 
introduce a diagonal density matrix p^ ag = xi <8 1 + (1 — 
X ) [ra <8> 1 + 1 <8> Pb] , where p A = R c \ +) (+| + R D \ -) (- 1 
and pb = Tc\+)(+\ + Td\ — )(— |, and a density ma- 
tri x p mt = R C T D \ +-)(- + |+ R D T C \ - +)(+'- I - 
~f^T c RcT D R D [e^ I -+>(- + I + e-* | +-)<+- 1] result- 
ing from the two-particle interference, with (f> a scattering 
phase. Here ® is a direct product between single particle 
density matrices at A and B and | — h) = | — )Ai\+)Bi 
with (-1 — I = (| — h))' etc. The orbital states \+}Ai/Bi 
(|— )Ai/Bi) denote the upper (lower) incoming leads to- 
wards detector regions A/B (see Fig. [3]). In agreement 
with Eq. decoherence 7 < 1, is introduced as a sup- 
pression of the two-particle interference |^r mt )(^r mt | — > 
p mt , where |$ int ) = VRc1d\ + -) - e^v^fel - +>• 
Following [lj we then write p P (E) = w p (E)o~ p where 
w p (E) = tr[ Pp (E)] = (1 - f) 2 fy[(R c T D + T C R D )(1 - 
X) 2 + 4x] and a p the normalized, energy independent 
density matrix of the emitted two-particle state. The 
entanglement production C P (E) = w p {E)C(a p ) is con- 
veniently quantified in terms of the concurrence C [llT |. 
ranging from for a separable state to 1 for a maximally 
entangled state. The total entanglement production in a 
time r, C p = (r/h) J dEC p {E), is then (Af = reV/h) 

C p = (AAi//2)max{47v / T Cj RcTr 1 i?z3-smh _2 x, 0}. (4) 

As shown in Fig. [2l C p decreases monotonically as a 
function of T, reaching zero at a critical temperature T£. 

Inserting the parameter values from the experiment, 
we find that C p w 0.LV and C(a p ) w 0.3, i.e. the state 
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FIG. 2: (color online) a) Entanglement production C p /N 
(blue, transparent) and C r /N (green, opaque) as functions of 
temperature kT/eV and coherence 7 for the semi-transparent 
2PI. b) Parameter Q as a function of kT/eV (blue line). Val- 
ues 0.25, 1CT 1 , 1(T 2 , 1(T 3 , 1(T 4 shown (grey lines). 



emitted by the 2PI is clearly entangled. However, this is 
a rather indirect deduction of the two-electron entangle- 
ment, an unambiguous demonstration of entanglement 
would be desirable. For this we turn to the reduced den- 
sity matrix p r , providing full information about the de- 
tectable two-particle correlations. 

We first express p r in terms of currents Iaoli^bp and 
low- frequency current cross correlators Q SacBP, with 
a, f3 = ±. Extending [5] to nonzero temperatures we find 



BP 



BP 



(Ve 2 /h) 2 2Ve 3 /h 



tr{[J& 



!bp\ 



Pr}- (5) 



The orbital current operators in the local basis {|±)}, 
including the rotations at the detector splitters, are 
J Aa = (1 + aS A (j z S\)/2 and I° p = (1 + pS B cr z S B )/2, 
with ay the Pauli matrix and Sa (Sb) the scattering ma- 
trix of the splitter at A (B). The results of Ref. [1(3] yield 



Pr = (1- H)[p A ® PB] +Hp v 



(6) 



Writing p r = uyay with w r = ix[p r ] — [RcTc + 
RjjTd](1 — H) + RcTjj + RdTc and ay the normalized 
density matrix we then define the total entanglement pro- 
duction during a time r as C r = Mw r C '(ay) ■ It is 

C r = 2J\f max{ ^T c RcT D R D [H(l + 7) - 1], 0}. (7) 

As C p , C r decreases monotonically with increasing T, 
reaching zero at T = T£, as shown in Fig. O 

Comparing Eqs. (T5]) and Q we find that C p > C r 
for Q(T) = ff/[4(l - H) sinh 2 x] < y/T c R c T D R D , in- 
dependent on 7 (see Fig. [2]). Consequently, for splitters 
away from the strongly asymmetrical, tunneling limit, the 
reduced entanglement constitutes a lower bound for the 
projected entanglement. In the tunneling limit, however, 
the reduced entanglement is larger than the projected 
one. This is the case since asymmetry has different ef- 
fects on C p and C r ; the critical temperature TJ, from 
Eq. ([7]), is independent on the splitter transparencies 
Tc,Td while T£, from Eq. ([J|, is reduced by increas- 
ing asymmetry. For the parameters in the experiment, 



Q(T) « 4 x 1CT 4 < ^R C T C R D T D « 0.25, showing the 
validity of the bound. However, C r w 0.01A/" and based 
on the measurement Q no conclusive statement can be 
made about C r and hence not about C p . 

A more detailed understanding of this finite tempera- 
ture readout problem can be obtained by comparing the 
properties of a v and oy. For perfect coherence 7 = 1 
and identical beam splitters Tc — To = T = 1 — 1Z 
one can (up to a local phase rotation) write o~ p / r = 
±£ p / r i® l+(l-£ P /r)|*s)(*s|, a Werner state [HI, with 
singlet weight Q 1 - £ p = 2KT sinh 2 x/[l + 2RT sinh 2 x] 
and 1 - Cr = H/(2 - H). Increasing kT/eV = 2/x from 
zero, ^p « 2e~ 2x /(RT) becomes exponentially small 
while Cr ~ 2/ai increases linearly. These qualitatively 
different behaviors, clearly illustrated in Fig. [IJ are a 
striking signature of how a small kT/eV, having negligi- 
ble effect on C(a p ), leads to a large suppression of C(ay). 

From Eqs. (H]) and ([7]) follows also a counter- intuitive 
result: finite amplitude of the AB-oscillations is no 
guarantee for finite two-particle entanglement. This 
is apparent for ay in the limit 7 = 1 and Tc = 
Td, since a separable Werner state, Cr > 2/3, can 
be decomposed [l3[ into a sum of product states as 
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atan[2/ A / s ] 



| + } +e «r[l-2n]/4 sin^/- !.) and6 » 
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?f /B = -acot[y A / B 
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a r d 
with 



^ = 

.A/B = 



{y/2^r + V3Cr = 2)/(V& ± V4 = 3£ r ) with +(-) for 
A(B). This classically correlated state gives, via Eq. ([5]), 
AB-oscillations with amplitude 2(1 — Cr)/(2 — Cr) = H. 

Moreover, the effect of decoherence, suppressing the 
two-particle interference, is similar for the projected and 
the reduced entanglement. In particular, for T = 0; 
Cp = C r — 2N^^/RcTcRdTd, finite for arbitrary strong 
dephasing 0] . This is a consequence of the 2PI-geometry, 
prohibiting scattering between upper (+) and lower (-) 
leads, i.e. pseudo-spin flips 0,0, 13|- Importantly, given 
the controllability of phase gates and beam splitters [IH , 
demonstrated in the electronic Mach-Zehndcr and 2PI 
experiments [1, [l6| , all the technics necessary for an en- 
tanglement test via reconstruction of p r are at hand. 

Another widely discussed [H, H, 0, [13] approach to de- 
tect the entanglement is to use a Bell inequality. Viola- 
tion of a CHSH-Bell inequality [l8| formulated in terms of 
currents and low-frequency current correlations demon- 
strates finite entanglement of p r . An optimal CHSH-Bell 
test demands the same number of measurement and level 
of experimental complexity as a tomographic reconstruc- 
tion of p r . From ay and ay, we can using [1 91 ] calculate the 
corresponding maximal Bell parameters SVf^p X and ft^ x , 
yielding for identical source splitters Tc = Td = T 



n™f r = 2 V / T+y (1 - Cp/r) 



(8) 



The CSHS-Bell inequalities are ft Bp/r < 2, with the lim- 
its = 2 for T = 1/2 plotted in Fig. [TJ It is clear 
that for the values kT/eV and 7 of the 2PI-experiment, 
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while Qsp < 2 in principle can be violated, a detection 
of entanglement by violating Qsr < 2 is not possible. 

Extending the above analysis to a general mesoscopic 
system, it can be shown [2(| that away from the strongly 
asymmetric limit (see below), C r > guarantees C p > 0, 
or cquivalently T£ > TJ, i.e. for a general conductor, 
nonzero entanglement of p r demonstrates that the emitted 
many-body state is entangled. This motivates a detailed 
investigation of the p r and C r in a general system. 

We consider a conductor, shown in Fig. [3l charac- 
terized by a scattering matrix 5* and connected via sin- 
gle mode leads to M > 2 reservoirs biased at eV and 
N — M > grounded. The conductor is also connected 
to detector reservoirs Act and B[3 via reflectionless split- 
ters with controllable phase gates. The splitter-phase 
gate structures perform local rotations, characterized by 
Sa and Sb, of the orbital states. All reservoirs are kept 
at temperature T. We assume linear response in ap- 
plied bias and S independent on energy E in the interval 
— kT < E < eV + kT of interest. The scattering is more- 
over assumed to be spin independent and we hence drop 
spin notation, presenting results for a single spin species. 




FIG. 3: (color online) Sketch of the general system. A meso- 
scopic scatterer is connected to M biased and N—M grounded 
source reservoirs and, via phase gates (j>A/B and beam split- 
ters, to grounded detector reservoirs A±,B±. The orbital 
basis states |±)xi/B< are displayed. For details see text. 

In order to treat general conductors we need to mod- 
ify Eq. © as S A aBp -> wn ere Sjf aB £ = 
S r J 4aB/3(eV, kT) — SAaBp(0, kT), the non-equilibrium ex- 
cess transport part of the correlator Q. Working in the 

orbital basis {| + +), | H — ), | — h), | )}, we get from 

scattering theory [8| the reduced density matrix 

Pr = t A t\ ® t B t B - H (t A t B <8> t B t\ ) P, (9) 

where P is a 4 x 4 permutation matrix with nonzero 
elements P tj = 1 for ij = {11,23,32,44}. The 2 x M 
transmission matrix t^ie) is the matrix amplitude for 
scattering from the M biased reservoirs to the two leads 
going out from the source towards region A (B). 

To obtain a compact expression for C r we make 
the singular value decompositions t& = UaTaVa and 
t B = UbTbVb, where the 2 x M matrix Ta = [ta,0] 
with ta = diag(-y/T^+, -\JTa-) and similar for Tb, and 



Ua,Va,Ub,Vb unitary. Inserting this decomposition 
into p r in Eq. ((9]) we arrive after some algebra at 



C r = 2M^/Ta+Ta-Tb + T b ~ max{F,0}, 
F = Hy/C^-y/il-HCMl-Hb) (10) 

where Cii C2 G [0, 1] are the eigenvalues of ZZ>, with Z a 
2x2 matrix with elements Z^ = (VlVs)ij for i,j = 1, 2. 

It can be shown *20] that for - - H)} < 

(T A l + T A X _ +T B l+T B 1 _)-\ C r > guarantees C p > 
0. This condition is expressed in terms of Ta±,Tb±, 
eigenvalues of the reduced single particle density matrices 
accessible by a reconstruction using average currents 0] . 

In conclusion, we have presented a theory for entangle- 
ment generation and detection in mesoscopic conductors 
at finite temperatures. It is found that under very gen- 
eral conditions finite reduced, detectable entanglement 
constitutes a witness for nonzero emitted, projected en- 
tanglement. The theory applied to the two-particle in- 
terferometer Q investigated experimentally by Neder et 
al 0] shows that while the emitted state is clearly entan- 
gled, the detectable entanglement is close to zero. Our 
results provide motivation for further experimental inves- 
tigations of entanglement in the 2PI. 
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